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We develop a finite-element technique that allows one to evaluate correction of the order of 
Gq to various transport characteristics of arbitrary nanostructures. Common examples of such 
corrections are weak localization effect on conductance and universal conductance fluctuations. Our 
approach, however, is not restricted to conductance only. It allows in the same manner to evaluate 
corrections to noise characteristics, superconducting properties, strongly non-equilibrium transport 
and transmission distribution. To enable such functionality, we consider Green functions of arbitrary 
matrix structure. We derive a finite-element technique from Cooperon and Diffuson ladders for 
these Green's functions. The derivation is supplemented with application examples. Those include 
transitions between ensembles and Aharonov-Bohm effect. 
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I. INTRODUCTION 

The theoretical predictions of weak localization |l| and 
universal conductance fluctuations |2| along with experi- 
mental discoveries in this direction 001 have laid a basis 
of modern understanding of quantum transport — trans- 
port in nanostructures — and have stimulated a consid- 
erable interest to the topic. Early studies mostly con- 
centrated on diffusive electron transport. Both effects 
arise from quantum interference that is described in the 
language of slow modes: Cooperons and Diffusons.0, Q 
Each mode of this kind brings a quantum (fluctuating) 
correction of the order of Gq = e 2 /irh to the classical 
Drude conductance G of the sample. This universal value 
sets an important division between classical conductors 
(G ^> Gq) where interference effects are small and quan- 
tum ones (G ~ Gq) where the transport is essentially 
quantum. 

A complementary approach to GQ-corrections comes 
from Random Matrix Theory(RMT) of scattering 0, 
This approach relates statistical properties of scattering 
matrix of a nanostructure to those of a certain ensem- 
ble of random matrices. GQ-corrections are understood 
in terms of fluctuations and rigidity of spectral density 
of these matrices. Although RMT approach can deal 
with diffusive systems, the most comprehensive setup in- 
cludes the so-called quantum cavity — an element whose 
scattering matrix is presented by a completely random 
unitary matrix of a certain ensemble. The cavity can be 
seen as a region of space where electron motion is suf- 
ficiently chaotic (either ballistic or diffusive) and where 
electrons can get in and out through some constrictions 
9j. The transport is determined by the propagation in 
the constrictions while random unitary matrix represent- 
ing the cavity is responsible for "randomization" of the 
scattering. The RMT approach does not necessary con- 
centrate on the total conductance. One can work with 
transmission distribution: the averaged density of eigen- 
values of transmission matrix squared. This transmis- 
sion distribution appears to be useful in a much broader 
physical context: it determines not only conductance of 



nanostructures, but also noise, full counting statistics of 
charge transfers, and properties of the same nanostruc- 
ture with superconducting leads attached @, 0] . It is a 
modern paradigm of quantum transport that an individ- 
ual nanostructure is completely characterized by a set of 
transmission eigenvalues while transmission distribution 
describes the averaged properties of random nanostruc- 
tures of the same design. This makes it relevant to study 
Gq corrections and fluctuations of transmission distri- 
bution [hH Il5| . The density of transmission eigenvalues 
is of the order of G/Gq, and Gq corrections are of the 
order of SG/G Q ~ 1. 

The microscopic Cooperon-Diffuson description is 
equivalent to a proper RMT approach. This is best illus- 
trated in the framework of more general supersymmetric 
theory \H\ that allows for non-perturbative treatment 
of fluctuations in quantum scattering. The Cooperons 
and Diffusons in this theory are fluctuations of super- 
symmetric field around the saddle point. For quantum 
cavity, only a single mode of these fluctuations is rele- 
vant. The inte grat ion over these modes reproduces RMT 
results [H [Hill. 

One can describe nanostructures in the framework of 
a simple finite-element approach usually termed "circuit 
theory" . The circuit theory has originated from attempts 
to find simple solutions of Usadel equations in super- 
conducting heterostructures [l6|. However, it has been 
quickly understood that theories of the same structure 
can be useful in much broader context: one can com- 
pute transmission distribution 0, , noise and count- 
ing statistics 0, investigate spin effects and non- 
equilibrium phenomena 20] . In circuit theory, a nanos- 
tructure is presented in a language similar to that of 
electric circuits: It consists of nodes, reservoirs and con- 
nectors. A node is in fact a quantum cavity, a connector 
can be of very different types: tunnel junction, ballistic 
contact and diffusive wire and is generally characterized 
by a set of transmission coefficients. In circuit theory, 
each node is described by a matrix related to the elec- 
tron Green's function. In the limit G ^> Gq the circuit 
theory provides a set of algebraic equations that allow 
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one to express the matrices in the nodes in terms of fixed 
matrices in the reservoirs. 

In this article, we present a technique to evaluate Gq- 
corrections for arbitrary nanostructures described by cir- 
cuit theory. To allow for various applications, we consider 
Gq corrections to multi-component Green's functions of 
arbitrary matrix structure. We are able to present Gq 
corrections in a form of a determinant made of deriva- 
tives of Green's functions in the nodes with respect to 
self-energy parts. In terms of a supersymmetric sigma- 
model, this corresponds to expansion of the action up 
to quadratic terms. However, the formulation we present 
does not contain any anti-commuting variables that com- 
plicate applications of sigma-models. The determinant is 
just that of a finite matrix, this facilitates the computa- 
tion of Gq corrections for nanostructures of complicated 
design. 

The structure of the article is as follows. To make 
it self-contained, we start with a short outline of cir- 
cuit theory of multi-component Green's functions ad- 
justed for the purposes of further derivations. In Sec- 
tion we derive microscopic expressions for Gq correc- 
tions and specify to finite-elements in Section ITvl Next 
Section [V] is devoted to description of spin-orbit scat- 
tering and magnetic-field de-coherence that are used to 
describe transitions between different RMT ensembles. 
Since Aharonov-Bohm effect plays an important role in 
experimental observation of Gq corrections, we explain 
how to incorporate it to our scheme in a short sepa- 
rate section IVII We illustrate the technique with sev- 
eral examples fSection lVIIf) concentrating on a simplest 
2x2 matrix structure that suits to calculate Gq cor- 
rections to transmission distribution. The examples also 
involve simplest circuits: One with a single node and 
two arbitrary connectors, a chain of tunnel junctions, and 
two-node four-junction circuit to demonstrate Aharonov- 
Bohm effect. We summarize in Section EnH 



II. FROM GREEN'S FUNCTIONS TO FINITE 
ELEMENTS 

In this article, we consider Green's functions of arbi- 
trary additional matrix structure with N c h indices. We 
do this for the sake of generality: This allows for descrip- 
tion of super-conductivity, incorporating Keldysh formal- 
ism and treating non-equilibrium and time-dependent 
problems. This is also extremely convenient, since the 
most relations in use do not depend on the "physical 
meaning" of the structure. We use "hat" symbol for op- 
erators in coordinate space and "check" symbol for ma- 
trices in additional indices. The Green's function thus 
reads: G = G(x,x'), x stands for (three-dimensional) 
coordinate. The general Green's function is defined as 
the solution of the following equation 



All physical quantities of interest can be in principle 
calculated from Green's functions. We address here the 
quantum transport of electrons in disordered media. In 
this case, one can work with common Hamiltonian 



H = H + u{x), H 



1 
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Vl + U(x). (2) 



t(x)-H G(x,x') = S(x-x'). 



(1) 



Here U (x) describes the design of nanostructure: Po- 
tential "walls" that determine its shape, and form bal- 
listic quantum point contacts, potential barriers in tun- 
nel junctions, etc. The potential u(x) is random: it de- 
scribes the random impurity potential responsible for dif- 
fusive motion of electrons, isotropization of electron dis- 
tribution function and, most importantly for this article, 
fluctuations of transport properties of the nanostructure. 
The physics at space scales exceeding the isotropization 
length (which is the mean free path in the case of diffusive 
transport) does not depend on a concrete model of ran- 
domness of this potential. Most convenient and widely 
used model assumes normal distribution of u(x) charac- 
terized by the correlator (u(x)u(x')) — wS(x — x'). It is 
important for us that both Hq and u(x) are diagonal in 
check indices. 

Most evident choice of self-energy matrix is T>(x) = — e, 
e being energy parameter of the Green's function. £ is 
more complicated in the theory of superconductivity. We 
will find it convenient, at least for derivation purposes, 
to work with arbitrary T,(x). We can also consider more 
general situation with non-local Yi(x,x'). 

Provided the conductance of the nanostructure is suffi- 
ciently high (G » Gq), one can disregard quantum Gq- 
corrections and work with semiclassical averaged Green's 
functions. Closed equations for those are obtained with 
the non-crossing approximation |2l|. They include the 
impurity self-energy, 

-t(x) - t imp - H G(e;x,x') = 5(x - x'), (3) 

t, imp (x) = wG(x,x) (4) 

At space scales exceeding the mean free path, one can 
write closed equations for the Green's function in the 
coinciding points [23, ■ It is convenient to change no- 
tations introducing dimensionless G(x) — > G(x,x)/iirv, 
with v the density of states at Fermi energy. For purely 
diffusive transport, one obtains Usadel equation 

V- j-mv[t{x),G{x)] = 0, 3 = -nD(x)isGVG, (5) 

D being the diffusion coefficient. The solutions of Us- 
adel equations are defined only if one takes into account 
boundary conditions at "infinity": Equilibrium Green's 
functions in the macroscopic leads adjacent to the nanos- 
tructure. It turns out that G satisfies unity condition 
G 2 = 1. For situations where the transport is not entirely 
diffusive, one would supplement the Usadel equation with 
boundary conditions of various kinds (c.f. 24]). 

An alternative way to proceed is to notice that the 
Usadel equation is almost a conservation law for the ma- 
trix current j , which allows for a finite clement approach. 
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This conservation law is exact at spa ce sc ale smaller than 
the coherence length estimated as V DH. It does not re- 
lay on assumption of diffusivity: It occurs because the 
Hamiltonian H commutes with the "check" structure. 
It is then natural to proceed to the separation of the 
nanostructure into regions where G(x) can be assumed 
constant. The next steps are the same as in traditional 
circuit theory of electric circuits, which exploits conser- 
vation of electric current and finite discretization of the 
system onto regions of approximately constant voltage 
— nodes. Each node is connected by connectors to other 
nodes or reservoirs, those representing macroscopic leads. 

After the discretization of the nanostructure we can 
write a Kirchhoff-like equation 



0. 



(6) 



where the indices c and a label the nodes and the con- 
nectors of the nanostructure. The current through a con- 
nector c which connects nodes cl and c2 reads 



Jr= T 



[Gel 7 G 



c2 



4 4^ 1 + T„ c (G cl G c2 + G c2 G cl - 2)/4' 



(7) 



{T°} being the set of transmission eigenvalues of the 
transmission matrix squared relative to the connector. 
Eq. [5] shows that the current j is not fully conserved; to 
this aim we included in Eq. a leakage current 

h,a = — 7— [S Q , G a ] , 5 a — ——, (8) 

where V a is the volume of the node; S a can be easily 
recognized as the average level spacing in the node. 

It is opportune to notice that the conservation law © 
can be obtained by requiring that the values assumed 
by the matrix Green's function in the nodes are such to 
minimize the following the action 



1 + — (G C \G C 2 + G c2 G c i — 2) 



(10) 



S Q = -^-TrT, a G a . (11) 

Oa 

The minimization of the action must be carried out pro- 
vided that the Green's functions in each node satisfy 
the normalization condition G^ = 1; this implies that 
the variation of the Green's function 8G a has to anti- 
commute with G a itself. A possible way to satisfy this 
condition is to write the variation as 



5G a = G a Sv a — Sv a Ga: 



(12) 



where no restriction is imposed on Sv a - The Kirchhoff 
equations are then re-written as 



£ = [<5,-S-] = 0. 

OV a SG a 



(13) 



The same action can be used at microscopic level too, 
even before averaging over u(x). To do this, we note 
once again that £(2:) form , at least formally, a set of 
parameters of our model. This can be straightforwardly 
extended to non-local operators Y,{x,x'). To this end, 
we define the action in terms of the following variational 
formula 



SS = - J dxdx'Ti [G(x', x)St(x, x')] 



(14) 



which is traditional in Green's function applications [2 ll 
l2fij . With making use of formal operator traces, this 
action can be written in terms of either £ or G, 



S = TV In £ + H 



S = Tt -In -G -(£ + JT)G 



(15) 
(16) 



apart from a constant. Eq. \I\is reproduced by varying 
either (|15J) or (|16(l over £ or G respectively. Since we are 
not planning to treat Gq corrections beyond perturba- 
tion theory, we do not attempt to use the exponent of 
this action as an integrand in some path integral repre- 
sentation of a result of exact averaging over u(x). This 
has been done in for supersymmetric matrices Q and 
in [2^| for Keldysh-based two-energy matrices Q(x . e, e'). 
We note, however, that the action used in [ll], |2a| is 
equivalent to if one substitutes G = Q. 

To this end, for averaged Green's functions the action 
in use is defined simply as 



SS = - 



J dxdx'Tr [< G(x',x) > St(x,x')] 



(17) 



III. Gq-CORRECTIONS to 
MULTI-COMPONENT GREEN'S FUNCTIONS 

In this Section, we outline a microscopic approach 
to Gq corrections suitable for multi-component Green's 
functions. The main idea is the same as in [THll^ where 
similar derivation has been done for purely diffusive case 
and for a concrete 2x2 matrix structure. It is known 
[H that diffuson and cooperon modes responsible for 
Gq corrections are presented as ladder diagrams made 
from averaged Green's functions. In usual technique, 
such diagrams have vertices corresponding to two cur- 
rent operators in Kubo formula. One does not have to 
work with vertices: Instead, one considers Cooperon and 
Diffuson contributions to the action that do not have 
any. The part of the action that presents Gq corrections 
is given by wheel diagrams. These wheels are made of 
either Cooperon or Diffuson ladder sections in a straight- 
forward way (see Fig. ^) 

We argue that the optimal way to present this part 
of the action is to double the existing " check" structure 
and to consider Eq. ^ in so-extended setup. Indeed, 
suppose we would like to address the most general ap- 
plication of Gq corrections: Parametric correlations [28j. 
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FIG. 1: Wheel diagrams that determine G^-corrections to 
the action. Single (double) line represents the Green's func- 
tions from the black (white ) block while dashed lines represent 
the averaging over disorder. The Diffuson (left) and Cooperon 
(right) wheel differ by mutual orientation of the lines. 



In this case, we start with two different "worlds" cor- 
responding to two different sets of parameters. In our 
case, all parameters can be in principle incorporated into 
S(:r). So we have "white" and "black" sets, T, b ,T, w , 
those define corresponding non-averaged Green's func- 
tions G bl G w . The averaging over the random poten- 
tial u(x) provides correlations between the two "worlds" 
and gives the rungs of the ladder diagrams that involve 
"white" and "black" Green's functions (Fig. [TJ. This 
is to obtain the Diffuson ladder wheel. The Cooperon 
ladder wheel is obtained by inverting the direction of 
the Green's function in one of the "worlds". This cor- 
responds to using transposed self-energies for this world, 
£m — * (E W ) T . This guarantees that the corresponding 
Green's functions are also transposed. Fluctuations at 
the same values of the parameters are naturally given by 
diagrams where and T* w either are the same in both 
"worlds" (Diffusons), or mutually transposed (Cooper- 
ons). 

Finally, we note that the doubled structure is also 
useful for evaluation of the weak localization correction. 
In this case, the last section of the Cooperon ladder is 
twisted before closing the wheel. 

To proceed further, let us introduce the operators K 
presenting a section of a corresponding ladder, 



b.cd 



K a 



(x,x') 



r>ab,cd t l\ 
cooper V ' ^ } 



w(x){G?{x,x')){G${x',x)) 
w{x)(G a b c {x,x'))(G b *{x,x')) 



(18) 
(19) 



where Latin letters represent "check" indices. As we have 
noted, the "white" Green's function is transposed for the 
Cooperon. Those are operators in the space spanned by 
the coordinates and the two check indices. 

Summing up all diagrams, we find the formal operator 
expressions for contributions to the action. For fluctua- 
tions, we have 



Tr 



ln(l - i\coopcr) 



Tr 



ln(l - K diB ) (20) 



diffuson and cooperon contributions are given by K^iff, 
K cooper respectively. For weak localization correction, 
one has to account for the fact that the last ladder section 



is twisted. We do this by introducing the permutation 
operator P, which exchange "check" indices, 



p j^ab.cd j^-ba,cd 



The contribution to the action becomes 



-Tr 



Pln(l -K, 



cooper/ 



(21) 



(22) 



The factor one half is included in the last formula to take 
into account that "black" and "white" Green's functions, 
in the case of weak localization, are not anymore inde- 
pendent. We note that K for the Cooperon is symmetric 
with respect to index exchange, so that K and P com- 
mute. The eigenfunctions and the eigenvalues of K are 
therefore either symmetric (K + ) or anti-symmetric (K~) 
with respect to permutations. We can rewrite the last ex- 
pression as a sum over these eigenvalues 



= \T. [Mi - K+) - Hi - k-)] 



(23) 



It is clear from the previous discussion that, in order 
to calculate the Gq corrections, one has to evaluate the 
eigenvalues of the ladder section K, both for Cooperon 
and Diffuson. We introduce now a method to compute 
this matrix easily. The observation is that the ladders 
under consideration are not specific for GQ-corrections: 
The same ladders determine a response of semiclassical 
Green's functions upon variation of £ pflj . 

To see this, let us go back to non-averaged Green's 
functions. We keep in mind that we have doubled the 
"check" space to include white and black sector. We add 
by hand a source term: The self-energy which mixes up 
"black" and "white" Green's functions, 5Y, bw (x). This 
source term will give rise to a correction to the Green's 
function in the same black- white sector. In the first or- 
der, we have 

5G bw {x,x') = 
- J dx 1 dx 2 G b (x } x 1 )S£ bw (x 1 ,x 2 )G w (x 1 ,x'), (24) 

which is best illustrated by the diagram in FigureEl Next 
step is to include the effect of the random potential u(x). 
We average Eq. |^| limiting ourselves to the non-crossing 
approximation and obtaining a set of ladder diagrams 
(Fig. |21) By summing up all the contributions we ob- 
tain the correction — taken in coinciding points — to the 
Green's function 



(SG bw (x,x)) = 



1 



^diff 



W{x) 1 _ k dii 



-St bw (x). (25) 



Eq. (|25|l is very valuable: it demonstrates that the re- 
sponse of the Green's function to the source term ST, bw 
is determined by the same ladder operator K, which we 
need to compute Gq corrections. 
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(5G bw ) 





■ 




FIG. 2: The response of Green's function SGbm on self- 
energy T,bw in semiclassical non-crossing approximation is de- 
termined by the sum of ladder diagrams. This allows to ex- 
press the wheel diagrams in terms of eigenvalues of response 
kernels. 



At the space scale of isotropization length K ~ 1. Usu- 
ally one is interested in the contribution arising from the 
larger space scale where Cooperon-Diffuson approxima- 
tion is valid. At this scale, the eigenvalues of K are either 
zero or very close to one. To see this, we cite the results 
for the homogeneous case with £ = const(a:). A conve- 
nient basis in "check" space is one where £ is diagonal, 
the eigenvalues being £„. The Green's function is di- 
agonal in this basis as well, G n = s a = sign Im(S a ). 
Owing to homogeneity, the section operator is diago- 
nal in the wave vector representation, its eigenvalues 
being K nm (q). Direct calculation similar to 6] gives 
K nm (q) = if s n = s m . lfs n ^s m , 



1 - K ab {q) » r Us b {^ a - £ b 



Dq 2 



(26) 



for Sr, ql <C 1(t = 2-kwv is the isotropization time). 
This equation makes the relation between our technique 
and the common technique for Cooperons and Diffusons 
in homogeneous media. Usually, the self-energy £ has 
equal number of eigenvalues with positive and negative 
imaginary parts. In this case, at each q, K has N 2 h /2 
zero and N 2 h /2 non-zero eigenvalues. 

Now we note that zero eigenvalues contribute neither 
to i|2(J|) nor to l|25[l . As to those close to 1, we may replace 
K by 1 in the numerator of l|25[l . We also note that SG 
can be presented as the derivative of the action (c.f. Eq. 
II 41) . Therefore, we can write the Gq corrections due to 
the diffuson modes to the action in terms of a determi- 
nant made of derivatives of the semiclassical action with 
respect to £ bu , , T, wb 



^G Q ,diff = - bidet' -w(x 



S 2 S 



(27) 



The " prime" sign of the determinant signals that the zero 
eigenvalues shall be excluded: det' is defined as the prod- 
uct of all non-zero eigenvalues. Indeed, as we have seen, 
some variations of self-energies do not change the Green's 
functions giving rise to zero eigenvalues. We also note 
that the Gq corrections are not affected by the concrete 
form of w{x): Since the determinant of matrix product 
is a product of their determinants, this matrix gives a 
constant contribution to the action which does not affect 
any physical quantities. 



IV. METHOD 

Let us now adapt the microscopic relation (|27|) to the 
finite-clement approach outlined in the Section [n] 

With all previous derivations, this step is easy. We just 
replace the actual cc-dependent G and £ by constants in 
each node. To get the action in these terms, one inte- 
grates over the volume of each node so that the formula 
(H1J) reads 



(28) 



where the summation is over the nodes. The discrete 
analog of the determinant relation l|27|l is now 



Sg q , diff = - In det' 



S 2 S 



V a §Y, wb 5'Ebw 



dndet' 



S 2 S 



const— In det' 



5 2 S 



(29) 



where we have introduced dimensionless response matrix 
f) a = iirY,/8 a and noticed that the matrix oc w brings a 
constant contribution to the action. The response matrix 
is determined from the solution of the Kirchhoff equa- 
tions at vanishing source term rf^ . It has N no d es • N 2 h /2 
non-zero eigenvalues and the same number of zero ones. 
We observe that at T, Wib — the eigenvalues of this ma- 
trix do not depend on the volume of the nodes, they are 
determined by the transmission eigenvalues of the con- 
nectors only, and are of the order of GJ Gq . Since rescal- 
ing of all conductances gives only an irrelevant constant 
contribution to the action, the Gq corrections depend 
only on the ratios of conductances of the connectors: This 
manifests the universality of these corrections. 

The circuit theory action Qis given in terms of G a . 
It is advantageous to present the answer for Sq q in 
terms of the expansion coefficients of the action around 
the saddle-point: The solution of semiclassical circuit- 
theory equations. That is, to use S 2 S /SG wb SG bw in- 
stead 5 2 S / 5f] wb 5fj bw . If the latter matrix were invertible, 
we would make use of the fact that S 2 S/5G wb 5G bw = 
(S 2 S / 8ij wb Sf) b . w )~ 1 . In fact, owing to the constrain G 2 = 
1, there is a big number of zero eigenvalues in the re- 
sponse matrix. So the task in hand is not completely 
trivial. 

We proceed as follows. We expand the action by re- 
placing each G in each node by 



G = Go + g - G g 2 /2 



(30) 



and collecting terms of the second order in g, fj. The 
form ll.''>l)|) satisfies the constrain G 2 = 1 up to second 
order terms provided gGg + Gog = 0. Let us work in 
dimensional space indexed with the "bar" 



the N 2 h N node 



index a composed of two " check" indices and one node in- 
dex, a = (a, b, a). We present the result of the expansion 
as 

5S = gfM- aW %»-r, & g%». (31) 

The variation of (fTTTfl under the constrain gGo + Gog = 
gives the response matrix 5 2 S/Sr]aSria- Next we consider 
matrix defined trough the following relation: 

- (g - G gG ) -> 

\ (g bw - G b g bw G w ) , (32) 

the last equation makes white-black block separation ex- 
plicit. We note that H a i is a projector. It separates "bar" 
space on two subspaces where g either commutes or anti- 
commutes with Go, and projects an arbitrary g onto anti- 
commuting subspace. Applying this projector to Eq. 1311 
we show that the projected matrix n^A/j^n^ is an in- 
verse of the response matrix within the anti-commuting 
subspace. 

S Gq =lndet'(n a5 M 5e n g5 )) = 

In det ;1U.\/,.1U + S- ah - U- ah ) (33) 

In the last equation we add the matrix 1 — II. This pro- 
cedure replaces all zero eigenvalues with 1, so one can 
evaluate a usual determinant. 

We remind that as far as fluctuations are concerned, 
there are two contributions of this kind coming from Dif- 
fuson and Cooperon ladder respectively. The weak lo- 
calization correction involves permutation operator that 
sorts out eigenvalues involved according to (|23|) . With 
this, the equations (|33|l and l|29|l give the Gq corrections 
in an arbitrary circuit-theory setup in the most general 
form. 



V. DECOHERENCE AND ENSEMBLES 

Until now we have assumed that the Hamiltonian com- 
mutes with the "check" structure and is invariant with 
respect to time reversal. This implies strict coherence 
of waves with different "check" index which propagate 
in the disordered media described by this Hamiltonian. 
Even small "check" -dependent perturbations of the sym- 
metric Hamiltonian give accumulating phase-shifts to 
these waves and may significantly change their interfer- 
ence patterns at long distances. Due to its random na- 
ture, such phase-shifts can be regarded as decoherence 
although this should not be confused with a real decoher- 
ence coming from interaction-driven inelastic processes 

In real experimental situations, two sources of such 
decoherence are usually of importance: spin-orbit scat- 
tering and magnetic field. Already early studies of Gq 
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FIG. 3: Second-order diagrams in spin-orbit interaction (left) 
and magnetic field (right) provide decoherence terms in the 
action. These terms describe transitions between pure RMT 
ensembles. 

corrections P, 3 have revealed their significant depen- 
dence on these two factors in the regime where those 
are too weak to affect the semiclassical transport. From 
the RMT point of view, these factors, upon increas- 
ing their strength, provide transitions form orthogonal 
ensemble (symmetric Hamiltonian) to two different en- 
sembles: symplectic (spin-orbit interaction) and unitary 
(magnetic field) [Tj. 

We show in this Section how to incorporate spin-orbit 
scattering and magnetic field into our scheme. The most 
convenient way is to present them as perturbative cor- 
rections to G-dependent action (Fig. |3) 

The spin-orbit scattering enters the Hamiltonian in the 
form H so = a a H a (x,x'), H a (x,x') = -H a (x',x), a a 
representing spin Pauli matrices in "check" space. In the 
second order in H a the averaging gives (Fig. |3J) 

S so = f dz-^-Tr [G{x)a a G(x)a a ] (34) 
J St so (x) 

At the level of microscopic approach, the spin-orbit scat- 
tering takes place anywhere in the nanostructure. In 
the finite-element approach, it is advantageous to ascribe 
spin-orbit scattering to nodes rather than to connectors. 
This is consistent with the main idea of our scheme: Ran- 
dom phase-shifts take place in the nodes. The spin-orbit 
contribution in each node a is obtained by integrating 
(|34Jl over the node, 

Sso = ^jpTr [G a a a G a &a\ (35) 

where rj so = ir / '(2tq So) 'S a ). 

The magnetic field is incorporated into the Hamilto- 
nian through the modification of the derivative, 

iV -» iV - -^rt, H A, 

ch 

where A is the vector potential and (T, 2 H = 1) de- 
scribes the interaction of different "check" waves with 
the magnetic field. In its simplest form, E# is the ma- 
trix in the white-black structure introduced such that 
f] b H = 1, — — 1 provided we describe a Cooperon. 
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This is consistent with the requirement that one of the 
Hamiltonians must be transposed to describe a Cooperon 
ladder. This is not the only plausible form of this ma- 
trix. For instance, in non-equilibrium superconductivity 
Ejy involves electron-hole Nambu structure. 

The magnetic field decoherence contribution can also 
be assigned to a node and reads 

S H = ^Tr [G a t H G a t H ] (36) 

where rjn = 7r/(r//5 Q ) and th corresponds to Cooperon 
magnetic decoherence time in a common theory. The 
latter is known to depend on the geometry of the node 
and its characteristics [3l]]. If the transport within the 
node is diffusive, 1/t h = 4(e / h) 2 D (A 2 ) where () denote 
averaging over the volume of the node. The vector po- 
tential is taken in the gauge where it is orthogonal to 
the boundaries of the node. We have order of magni- 
tude, 1/t h 6 ~ ($/$o) 2 (G„ ode /G Q ) $ being magnetic 
flux through the node, $o = Trh/e being flux quantum, 
and G no de being a typical conductance of the node. The 
latter is limited by its Sharvin value in the ballistic regime 
where the isotropization length is of the order of the node 
size. 

The magnetic field produces not only random but 
also deterministic phase-shifts. This gives rise to the 
Aharonov-Bohm effect discussed in the next Section. 

To find the effect of decoherence terms (|35fl and l|36[) 
on the eigenvalues forming the localization correction, 
we expand the action as done to obtain Eq. 1311 The 
decoherence contribution to M is diagonal in node in- 
dex, and can be made diagonal in "bar" index by proper 
choice of the basis in "check" space. For instance, if 
no external spin polarization is present in the structure, 
spin-orbit contribution is diagonal in the basis made of 
singlets and triplets in spin space. The simple realization 
of Efl- mentioned is automatically diagonal. If in addi- 
tion this diagonal contribution is the same in all nodes, 
both decoherence effects just shift the eigenvalues of M 
corresponding to the symmetric Hamiltonian. This gives 
an extremely convenient model of decoherece effects. 

The action for fluctuations is modified as follows: 

Smfj = Yj ln ( M «) + 31n ( M « + Vbo) (37) 

n 

Scooper = 22 ln(M„ + i] H ) + 3 ln(M n + i] so + r) H ) (38) 

n 

where summation goes over non-degenerate eigenvalues 
of M and factor 3 comes from three-fold degeneracy of the 
triplet. To derive the modification for the week localiza- 
tion contribution, we note that singlet and triplet are re- 
spectively antisymmetric and symmetric with respect to 
permutation. Therefore, triplet extensions of symmetric 
and antisymmetric eigenvalues are respectively symmet- 
ric and antisymmetric. The weak localization correction 
thus reads 




FIG. 4: Ahronov-Bohm effect: the topological phase arising 
from the magnetic flux $ is incorporated into the boundary 
conditions that relate the wave functions on two sides of an 
imaginary cut. The cut can be made in any place. 

AI + (~^ being (anti)symmetric eigenvalues of M. 

Since eigenvalues of M are of the order of G/Gq, 
the decoherence effects become important at rj so ,rjn — 
G/Gq, that is, when inverse decoherence times match 
Thouless energy E t h = (G/Gq)S of the node, 
1/t so , 1/t h ~ E th - 

VI. AHRONOV-BOHM EFFECT 

Aharonov-Bohm effect plays a crucial role in experi- 
mental observation and identification of Gq corrections 
(see e.g. ( 32]). Therefore it must be incorporated into 
our scheme, and in this Section we explicate how to do 
this. This extends the results of Ref. [33 where AB ef- 
fect was considered in superconducting circuit theory. In 
the following, we do not consider any orbital effects of 
the magnetic field but just the topological one. 

Let us suppose that the nanostructure presents a closed 
ring threaded by a magnetic flux <f>. As explained above, 
in presence of magnetic field the momentum operator has 
to be modified according to 

iV -> zV - e/ht H A, 

where A is the vector potential. Neglecting orbital ef- 
fects, one can get rid of the vector potential in the 
Schrodinger equation by a gauge transformation. Let us 
have a ideal cut in the nanostructure that breaks the loop 
(Fig. 0]) The topological effect of the flux can be incor- 
porated into a boundary condition for the wave-function 
^l,r on two sides of the cut, tpi = e~xjp(icj)AB^H)TpR- The 
phase of the wave function therefore presents a disconti- 
nuity at the cut that is equal to ±4'ab,4 > ab — n<f>/<&Q. 
Since the transformation does not explicitly depend on x, 
it can be immediately extended to semiclassical Green's 
functions. So that, those functions at two sides are re- 
lated by 

G L = exp(i4> AB t, H )G R exp(-i(j>AB'EH)- (40) 



This solves the problem at microscopic level. Once the 
nanostructure has been discretized to finite elements, we 
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note that the cut always occur between a connector and 
a node. The most convenient way to deal with the gauge 
transformation 1)40(1 is to put it into the action of the 
corresponding connector. To do this, we observe that 
the Green's function at the right end of the connector is 
not Gr of the node anymore: since the cut is crossed, it 
is eventually Gl given by H4U|) . The connector action in 
the presence of flux is therefore 



s„ 



T n - 

1 + — {G c iG C 2{4>ab) 



+G c2 {<Pab)G c1 - 2)] (41) 



where 



G{(j) A B) = exp(i(j) AB t H )Gexp(-i^ AB t H ) 



(42) 



One checks that the variation of the so-modified action 
reproduces the Kirchhoff laws for matrix current given in 
[33j. Owing to global gauge invariance, it does not matter 
to which connector and to which end of the connector the 
Aharonov-Bohm phase is ascribed. If the are more loops 
in the nanostructure, more connector actions have to be 
modified in such a way. 



VII. EXAMPLES 

In the previous sections, we operated with general ma- 
trix "check" structure to keep the discussion as general 
as possible. In this Section, we will give a set of examples 
to illustrate concrete applications of the technique devel- 
oped. For the sake of simplicity, we choose the simplest 
matrix structure that gives a sensible circuit theory. We 
consider 2x2 matrix Green's functions whose values in 
two terminals can be parametrized by a single parameter 







(43) 



The use of these matrix structures is that they give 
access to a fundamental quantity in quantum transport: 
Transmission distribution of transmission eigenvalues of 
two-terminal nanostructure. We have to explain this re- 
lation before going to concrete examples. The averaged 
transmission distribution is defined as 



p(T) = y(6(T-T n )) 



(44) 



where the sum is done over all the transport channels 
and the average is, in principle, to be intended over an 
ensemble of nanostructures of the same design. For G ^> 
Gq the transmission eigenvalues are dense in the interval 
[0, 1] and self-averaging takes place. 

Let us take a connector and set Green's functions at 
its ends to G(0i), G(</>2). By virtue of lfTU|l the connector 
action reads 



s{4>\ - 4>2 



1 - T„ sin 2 ■ 



(45) 



Ci 



FIG. 5: The chain of tunnels junctions of the same con- 
ductance Gt- In our finite-element approach, there are two 
Cooperon and two Diffuson eigenvalues per node. In the limit 
of N — > 00 the chain is shown to be identical to a continous 
quasi-one-dimensional diffusive wire. 



The trick is to regard the whole nanostructure as a 
single complex connector between left and right reser- 
voir, and set the Green's functions in the reservoirs to 
G(0),G(cj)). The total action © becomes now the con- 
nector action of the whole nanostructure and defines the 
transmission distribution in question, 



S(</>) 



dTp{T) In 



1 -Tsin^- 



(46) 



If one computes the ^-dependence of 5, the transmis- 
sion distribution can be extracted from its analytic con- 
tinuation on complex cj) ([lfl|) 



P(T) = 



1 



-.Re 



dS , 



2icosh~ 1 -= - 0+) 



T 



(47) 

Circuit theory of the Section [H] gives the answer in 
the limit G 3> Gq. The weak localization contribution 
S w i gives Gq correction to the transmission distribution. 
The fluctuation contribution Sg q — Sdi // + S CO oper that 
depends on two parameters c Wft, gives correlations of 
transmission distributions (c.f. 15]) 



J dTdT 1 < p{T)p{T') > 



hi 



1 — 1 sin — 



In 



1 — 2 sin — 
2 



(48) 



A simple application of the above formulas are Gq correc- 
tions to the conductance. Those are given by the deriva- 
tives of corresponding actions at <f>b,w = 0, 



SG w i _ d 2 S w i{4>) 
G Q ~ dcp 

= 4 



<C GbG w > „ d A S GQ (f>w) 



(49) 
(50) 



A. Junction Chain 

The first example is a chain of tunnel junctions. We 
will study Gq corrections for a chain of N identical junc- 
tions that connects two reservoirs (Fig. [5J) The connec- 
tor action for a tunnel junction assumes a very simple 
form: S c = (G T / '4G Q )Tr (G lc G 2c ) , G T being conduc- 
tance of the tunnel junction. Several tunnel junctions in 
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series, however, provide a good approximate for diffusive 
wire. Therefore, in the limit N — ► oo we can compare 
Gq corrections with the known results [T3 . ITHj l for cor- 
rections to transmission distribution of a one-dimensional 
diffusive conductor. 

We set Green's functions in the reservoirs on the left 
and on the right to G(0),G(<p) (c.f. Ea l43fl respectively. 
The semiclassical action for the system reads 



S 



„ N-l . N-l 

( ' r TrG fc G fe+1 - f ^ TrE k G k , 

fc=0 s fe=l 



4G 



(51) 



Here /c = 1, . . . , N — 1 labels the nodes, while k = and 
k = N identify left and right reservoir respectively. All 
the nodes are assumed to be identical with the same 5s- 
Since all the junctions are identical, the semiclassical 
solution is easy to find: the "phase" <j> drops by the same 
amount at each junction, and the solution reads Gk — 
G(kcj)/N), provided Efc oc 1 within each 2x2 block. This 
gives the optimal value of the action, 



_ NG T 

s= 2G^ C0 \n 



(52) 



from which one can evaluate the semiclassical transmis- 
sion distribution by using relation l|47[) . In the limit 
N -> oo S = -(G D /4:GQ)<t> 2 , G D ee G T /N being 
the conductance of the whole chain. This gives well- 
known transmission distribution for diffusive conductor, 
p(T) = (G d /2G q )/TVT^T [13. 

As explained above, to calculate Gq corrections we 
augment the "check" dimension of the Green's functions 
by introducing the "black" and "white" structure. Con- 
sequently, the parameter gets a "color" index b or w. 
The semiclassical solution for resulting 4x4 matrix is 
non-zero in bb and ww blocks, 



^ _ / G(kfa/N) 





&{k<j> w /N) 



(53) 



Now we shall derive the matrix M eigenvalues of which 
determine Gq corrections. It is advantageous to use a 
parametrization of the deviations from the semiclassical 
solution,^, which automatically satisfy gG + Gg — in 
each node. To this end, we rewrite the action l|51|l in a 
special basis: That one where G° are diagonal in each 
node, 



10 
0-100 
10 
-1 



(54) 



Then the deviation of the form 

n „bw 



gk 



( 
o 



V 9t% 



o 







/ 



(55) 



satisfies the condition mentioned. 

In this basis (see Appendix [3] for details) the action 
reads 



„ N-l . N-l 

S --^r- ^TrGfeMfe+iL- 1 -^ ^TrS fc G fc) (56) 



2G C 



fc=0 fe=l 

where bb (ww) block of L is given by 

T bb( ww ) _ ( cos(/W/A) ism(ct> b W/N) \ 

\ism(cl) b ^/N) cosi^/N) J' 1 ' 

We expand the Green's matrices according to Eq. (I3U|> . 
write the quadratic form in terms of g and diagonalize it 
(see Appendix^) to find the following set of eigenvalues, 
(1 = 1. ..N-l) 



4G^ 
Gt 



-M, 



±rj.b 



x , r r *i 

) = 2 cos — — cos —— cos — — cos — 



2N 



2N 



N 



■ » ■ 2 4> b ■ 2 4> W o^l 

cos — T 1/4 sin sin cos 2 he. 

2V V 27V 2N N 



N 



(58) 



e ee 27rGQ(E b — Y, w )/GTiSs measures the difference of 
Green's function energy parameter in bb and ww blocks in 
units of a single-node Thouless energy. To obtain eigen- 
values that determine the weak localization contribution, 



we set 



, e = 0. This yields 



M+.U) = —Ilcosf— ) 
wl ' lKV > AG Q y N J 



nl 
N 



Gt 
4Gc 



kI 
N 



N 



(59) 



(60) 



Let us discuss the weak localization correction first. If 
we neglect decoherence factors, we can sum up over / to 
find a compact analytical expression 



S w i = ln(cos— , 



1 



In 



2^ ' 



(61) 



In the limit N — > oo this reproduces the known correction 
for one-dimensional diffusive wire |l5j |. 



S w i = ^ In 



(62) 



It is interesting to note that the weak localization correc- 
tion is absent for N = 2. We will see below that this is 
a general property of a single-node tunnel-junction sys- 
tem. It was observed in [l5j that a part of weak local- 
ization correction in diffusive conductors is universal: It 
depends neither on the shape nor on the dimensionality 
of the conductor. The universal part is concentrated near 
transmissions close to one and is given by 



Spwi(T) 



1 



6(T-1), 



(63) 
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while the non-universal part is a smooth function of T. 
The relation (|MT)) possesses this property at any N, since 
the universal part comes from the divergency in 161(1 at 
(\> = 7r where the eigenvalue M~ t , goes to zero. Our ap- 
proach proves that this correction is universal for a large 
class of the nanostructures, not limited to diffusive ones, 
for any nanostructure where transmission eigenvalues ap- 
proach 1. This is guaranteed by the logarithmic form of 
the action. If M~ t 1 oc (ir — (j>) at cj) — > 7r, the correc- 
tion is given by ((63(1 irrespective of the proportionality 
coefficient. 

Expanding 1(61(1 at </> — > we find the correction to the 
conductance of the tunnel junction chain, 



5G U 



Gq 



1 (N 
3 



l)(iV-2) 



N 2 



(64) 



This is written for orthogonal ensemble, a well-known 
factor (1 — 2//3) defines the correction for other pure en- 
sembles. The effect of spin and magnetic decoherence 
can be taken into account by shifting the eigenvalues 1(591 
160(1 according to Eq. since the decoherence factors in 
each node are the same, (h,so = Vh,so(^Gq/Gt)- 

The S w i is still given by analytical although lengthy ex- 
pression (See Appendix A) . The correction to the trans- 
mission distribution corresponding to this expression is 
plotted in Fig. [fJlfor different strengths of spin-orbit cou- 
pling to illustrate the transition between the orthogonal 
and simplectic ensembles. The correction to the conduc- 
tance is given by 



SG 



wl 



Gq 



= 3F(NXh + (so)-F(N,(h) 



(65) 



where we define an auxiliary function F(x, N) 



F(x,N) = - 



(N-l) 1 {2 + xf 



■is 



N 2 N 2 x{A + x) 
1 2 2N + (2 + x + y/x(4 + x)) 

2N 



Let us discuss the parametric correlations. Without 
decoherence factors and at the same energy (e = 0) one 
can still sum up over the modes to obtain an analytical 
expression, 



Sniff = Scoop = (N - 1) In ( cos — + cos — 



sm 



■t>b±4 



(66) 



± \ k,iii 2N 

The fluctuation of conductance obtained with 1(50(1 reads 



((SG 2 )) 
r 2 



2 N 4 + 15N - 16 



15 



TV 4 



(67) 



and converges to the known expression for quasi-one- 
dimensional diffusive conductor at N — > 00. We notice 
that this convergence is rather quick, the fluctuation at 
N = 5 differs from asymptotic value by 10% only. We see 
thus that the diffusive wire, that in principle contains an 
infinite number of Cooperon and Diffuson modes, can be, 
with sufficient accuracy, described by the finite-clement 
technique even at low number of elements. 

Another point to discuss concerns the correlations of 
transmission eigenvalues T n , those can be obtained by an- 
alytic calculation of ((48(1 . It is instructive to concentrate 
on the relatively small eigenvalue separations those are 
much smaller than 1 but still exceed an average spacing 
~ Gq/G between the eigenvalues, Gq/G < \T-T'\ < 1. 
We observe that the correlation in this case is determined 
by the divergency of S at (j>b — 4> w — > ±2ir. Indeed, M-f 
approaches in this limit. This again suggest the univer- 
sality of these correlations. Indeed, as shown in 0] for 
diffusive conductors, the correlations in this parameter 
range are determined by universal Wigner-Dyson statis- 
tics and reduce to 



((p(T)p(T0)) = -^Re^l 



(68) 



Since the conductance fluctuations are contributed by 
correlations of T n at scale ~ 1 as well, they are not uni- 
versal. We plot in Fig. the correlator of conductance 
fluctuations as function of energy difference at several N. 



B. AB ring 

In this subsection we exemplify evaluation of AB effect 
within our scheme. We concentrate on the simple circuit 
presented in Fig. [S] It contains four tunnel junctions and 
two nodes labeled A and B. The conductances of the 
junctions are chosen to re-use the results of the previous 
section for a chain of three tunnel junctions: The solution 
of semiclassical circuit theory equations is given by lf5^|) 
for N = 3. The action reads 

G T f - - 1 - - 
S = —^Tr G L G A + ^GaGb 
4Uq ( 2 

+~G A G B {cj> AB ) + G B G R \ -i~ TrSjGj. (69) 

' i=A,B 

where Gl.r are Green's functions in the reservoirs and 
Gb{<I>ab) is modified according to Eq. |JJJ To study 
correlation of conductance fluctuations, we consider dif- 
ferent Green's functions for white and black block, and 
subjected to different fluxes 4>\ B 7^ F° r weak lo- 

calization correction, we set (f)\ B = —4>ab = ^ab 

To calculate the matrix M we use again the basis where 
G' ' are diagonal and the parametrization for g intro- 
duced in the previous subsection. It reads 



M 



M d M od 
M* nd M d 
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FIG. 6: Weak localization correction to the transmission dis- 
tribution of a system of four identical junctions at different 
values of spin-orbit parameter (so- We plot here cumulate 
correction X{T) = J^dT'T'8p{T'). X(l) represents the uni- 
versal singular part of the correction (c.f. Eq. 1631 while X(0) 
gives correction to the conductance. The lowest curve cor- 
responds to strictly zero (so and therefore represents pure 
orthogonal ensemble. Its negative value at T — 1 is partially 
compensated by positive non-universal contribution coming 
from T ~ 1 so that the resulting correction to conductance 
5G w i/Gq = X(0) ft! 0.2. The two higher curves correspond 
to relatively small values of (so, 0.05 and 0.4. While they are 
close to the orthogonal ensemble result at T ~ 1, their behav- 
ior at T ~ 1 is quite different: the universal correction is that 
of symplectic ensemble and is of positive sign. The highest 
curve corresponding to (so = 10 is close the cumulate correc- 
tion of pure symplectic ensemble, X sym (T) = — X ort (T) /2. 




R 



G T /2 

FIG. 8: A simplest circuit to exemplify AB effect consists of 
two nodes and for (tunnel) connectors. AB phase modifies 
the Green's function on the right side of the lowest connector 
with respect to Gb- The cut is given by dotted line. 



where 2x2 blocks Md,od are given 



G T M d 
4Go 



— cos ip — cos ip + e 




■ COS if 
2 





- cos <p u 



sin %r sin % 

2 A 



The parameter e which characterizes the energy differ- 
ence between black and white Green's function is defined 
as in the previous subsection. At e = we obtain explicit 
expression for the Diffuson eigenvalues (the Cooperon 



ones are obtained by 



and 



'AB 



J AB 



o 

to 



0.15 



0.1 



0.05 



4 5 6 

E/Et h 





N=5 




N=4 




N=3 




N=2 


* \ ■ 
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FIG. 7: Energy dependence of the correlator of conductance 
fluctuations <C G(0)G(E) ^> for chains with different number 
of junctions N. The energy difference is normalized to Thou- 
less energy of the whole chain, E t h = S s Gt /2-kGqN 2 . Note 
the fast convergence of the correlator to that of diffusive wire 
for big N and negative correlations at large E. 



4Gq 
Gt 



M 



± cos 



± cos 



i>AB 



1.2 



AB 



n 1,2 

(jr J 1 



COS 



i>AB 



,J AB 



<t>b 



() 



(70) 



(71) 



The weak localization correction to the action reads 

(72) 



1 ^cos 2 ^/3)(4-cos 2 (^b)) 



4 cos 2 (0/3) 



cos^ 



-'AB J 



from which we calculate the correction to conductance as 
function of the flux 



SG 



wl 



4 cos 



0AB 



Gq 



9(7-cos(20ab)) 



(73) 



We see that the weak localization correction cancels at 
half-integer flux 4>ab = it. This is because the junctions 
forming the loop are taken to be identical. The flux de- 
pendence exhibits higher harmonics indicating semiclas- 
sical orbits that encircle the flux more than once. 
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FIG. 9: The simplest possible circuit comprises one node, 
two connectors(L and R), and two reservoirs where Green's 
functions Gl r are fixed. 



For the correlator of conductance fluctuations we ob- 
tain 



(5G 2 



E 



259 - 4cos( 



AB 



± 



cos2( 



AB 



± 



81(cos( 



'AB 



± 



If 



(74) 

where plus and minus signs indicate Cooperon and Dif- 
fuson contributions respectively. The higher harmonics 
are present as well. 



C. Two connectors and one node 

Probably the simplest system to be considered by cir- 
cuit theory methods consists of a single node and two 
connectors (Fig. EJ). Since in this case there are only 
N c h eigenvalues, one can straightforwardly elaborate on 
complicated arbitrary connectors. For this setup we are 
still able to find an analytical expression for Cooperon 
and Diffusion eigenvalues. This allows us to get an ex- 
pression for the weak localization correction to the con- 



ductance which was vanishing in the case of two tunnel 
junctions. Each connector is in principle characterized by 
the distribution of transmission coefficients {T^}, {T^}, 
or, equivalently, by the functional form of the connector 
action given by Eq. QlUjl. Sl,Sr. ■ The action for the 
whole system reads 



S = S L (G L ,G) 



S R (G,G R )-i-Tr£G. 



(75) 



G being the Green's function of the node. We employ 
2x2 matrices parametrized by (|43l) , and set the Green's 
functions in the left and right reservoir to G(—(j)/2) and 
G((f>/2) respectively. The saddle-point value of G is given 
by the phase x an d for a general choice of Sl and S R does 
depend on 0, x = x{4>)- The total action in the saddle 
point is therefore S(<p) = S L (x + <P/2)+S R {x-<P/2). We 
expand the Green's function according to Eq. [3U] The 
second order correction to the action in this case reads 

S {2) = -i-^TrEa 
5 S 



2 v W 



9 2 G°G L 



+ ^({G°,G L }-2)] 



T 1 



1 2 



4[l + ^({G°,G L }-2)] 



(gG L gG L +g 2 )+(L^ R) 



Acting like in the previous subsections, we find two fol- 
lowing Diffuson eigenvalues (the Cooperon ones are ob- 
tained by the substitution <f> w — > —<fr w ). 



AI 



) = E 



cot[x(^ 



l-cos[ X (0 b ) + ^Tx(r)T^] 



1_cos[x(^)-£tx(0±£] 



E 

i—b,w 

E 



i—b,w 



1 



sm^[ X (^)-072] 



X'{<t> 1 ) + 1/2 V ' 1 V ; 2 
I'id) 1 ) 6 i 



r 



(77) 



Here we introduce I(<jj) = dS/dtfi to characterize the 
derivative of the total semiclassical action. We see that 
M~ approaches zero in the limit fibifiw ~> ±tt, pro- 
vided I{<f>) stays finite. As discussed, this divergency 
guarantees the universality of the correlations of trans- 
mission eigenvalues. 



1. Symmetric setup 

If we set S R — Sl, x(</> 6 ^) is zero regardless the 
concrete form Sl- The total action therefore reads 
S((f>) — 2Sl(4>/2). The eigenvalues (|77j> take a simpler 
form. To compute the weak localization correction to the 



action we set <b = 



to find 



Below we specify to three different cases. 



S w i = - In 



^tan^ 

m 2 



const. 



(78) 
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For tunnel junctions, I(<fi) oc sin(0/2) and the correction 
disappears. 

Expanding Eq. 1781 in Taylor series near (f> — * we 
reproduce the well-known result of 13] for weak local- 
ization correction to conductance. To comply with the 
notations used there, we characterize connectors with 



n n' 



5G 



W I 



G c 



h_ 

'4*i' 



(79) 



Similar expansion of Diffusion and Cooperon eigenvalues 
reproduces the result of ^| for conductance fluctuations, 



((SG) 2 ) = 3tl + 2tf-2ti(t 2 +*3) 
Gq ~ 8t? 



(80) 



2. Diffusive connectors 

It is instructive for understanding circuit theory of Gq 
corrections to specify the relation (|78|l to diffusive con- 
nectors. Since in this case I((f>) oc <fi, we obtain 



Swi,node = \; ln(tan(<£/2) /<f>) 



(81) 



Two connectors, single node situation can be easily real- 
ized in a quasi-one-dimensional wire with inhomogeneous 
resistivity distribution along the wire. A low-resistivity 
region would make a node if bounded by two shorter re- 
sistive regions that would make the connectors. On the 
other hand it has been proven in j2?J that the weak lo- 
calization correction in inhomogeneous wires does not de- 
pend on the resistivity distribution. Therefore, it has to 
be universally given by S w i t id = (1/2) ln(<j>/ sm(p) ^ 
S w i,node- How to understand this apparent discrepancy? 

This illustrates a very general point: Gq corrections 
may be accumulated at various space scales ranging from 
mean free path to sample size. The experimental obser- 
vation of the corrections relies on the ability to separate 
the contributions coming from different scales, e.g. by 
changing magnetic field 4] . With our approach, we eval- 
uate the part coming from interference at the scale of the 
node. The part coming from interference at shorter scale 
associated with the connectors is assumed to be included 
into transmission distribution of these connectors. 

For our particular setup, this extra contribution comes 
from two identical connectors. Since only half of the 
phase 4> drops at each connector, the contribution equals 
2iS u ,z j id(^/2). Summing up both contributions, we obtain 

Swl,node + ^S w l y id{4>/2) = S w l yt ot = S w l y id(4>) 

That is, the weak localization correction in this case re- 
mains universal provided the contribution of the node is 
augmented by contributions of two connectors. 



3. Non-ideal Quantum Point Contact 

The transmission distribution of an ideal multi-mode 
Quantum Point Contact (QPCj with conductance Gb ^> 
Gq is very degenerate since all T n = 1 or 0. This degener- 
acy is lifted if the QPC is adjacent to a disordered region, 
even if the scattering in this region is weak. This can be 
modeled as a connector with conductance Gd 3> Gb in 
series. 

The weak localization correction to the conductance 
has been calculated already a while ago ( [3{| ) . In the rel- 
evant limit, it is parametrically small in comparison with 
Gq, SG w i — —Gq(Gb/Gd)- Usual way to verify the ap- 
plicability of semiclassical approach to quantum trans- 
port is to compare the conductance of a nanostructure 
with the weak localization correction to it. For generic 
nanostructure, this gives G 3> Gq. However, for our par- 
ticular example SG w i <C Gb even for a few-channel QPC 
where Gb — Gq. So the question is: Is semiclassical 
approach really valid at Gb — Gq? 

To answer this question, we compute the weak local- 
ization correction to the transmission distribution. Since 
system is not symmetric, we make use of the full expres- 
sion (|77l) . In the limit of Gd 3> Gb, the relevant values 
of <j> are close to ir. We stress it by shifting the phase, 
fi = ix - (f>, \fi\ < 1. 

The circuit theory analysis in semiclassical limit gives 

M 



7T U „ 

*=2 + VT + i?c ' 



where R c = AGb/Gd <C 1. This gives to the following 
distribution of reflection coefficients (c.f. Eq. I47|l 



We use the above relations with (|77|l to find the 
Cooperon eigenvalues, 



M+{p) = -G 



D 



/i/4 \ \ + R c 



^ + R C 



Rc/2 



,(82) 



(83) 



This yields the weak localization correction to the cur- 
rent, 

2R C 



Iwl (//) 



fi{ f i 2 + AR C ) • 



The resulting correction to the transmission distribution 
consists of two delta-functional peaks of opposite sign, 
those are located at the edges of the semiclassical distri- 
bution, 

8 Pwl (JZ) = I [S(R -R c )- 8{R)] . (84) 
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To estimate the conditions of applicability, we smooth 
the correction at the scale of R c . This gives \Sp\/p ~ 
Gq/Gb and the semiclassical approach does not work 
at Gr ~ Gq. This agrees with RMT arguments given 
in |30j . The correction to the conductance calculated 
with (|84|) agrees with the value cited above and indeed 
is anomalously small. 



Ahronov-Bohm effect. We obtain general results for a 
single-node system with two arbitrary connectors and 
check their consistence with the well-known results for 
quantum cavity. This allows us to improve understand- 
ing of quantum interference in inhomogeneous diffusive 
wires and non-ideal quantum point-contacts. 



VIII. CONCLUSIONS 

We present a finite-element method to evaluate quan- 
tum corrections — typically of the order of Gq - 
to transport characteristics of arbitrary nanostructures. 
This includes universal conductance fluctuation and weak 
localization. We work with matrix Green's functions of 
arbitrary structure to treat a wider class of problems that 
includes superconductivity, full counting statistics and 
non-equlibrium transport. At microscopic level, the cor- 
rections are expressed in terms of Diffuson and Cooperon 
modes of continuous Green's functions. We employ a 
variational method based on an action to formulate a 
consistent finite-element approach. 

We illustrate the method with a set of simple and phys- 
ically interesting examples. All examples are based on 
2x2 matrices, this suffices to calculate the transmission 
distribution of two-terminal structures. We show how a 
chain of tunnel junctions approaches the diffusive wire 
upon increasing the number of junctions and study tran- 
sitions between ideal RMT ensembles in the chain. We 
consider the simplest finite-element system that exhibits 
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APPENDIX A: JUNCTION CHAIN 

In this Appendix we illustrate how to find the eigenval- 
ues of the matrix M defined in section llVl for the chain 
of tunnel junctions introduced. We consider the action 
in the rotated basis presented in text 



S 



Gt 
4Gc 



N 
i=0 



TrGiLGi+rL- 1 - ^ V Tr£<(?<. (Al) 



JV-l 



1=1 



We expand the Green's matrices according to Eq. (I3(jp 

The second order correction to the action reads (the 
first order is zero because we perform the expansion 
around the stationary point) 



J 



5(2) = E Tr {»■'■»■■ i '■ ' - l^LGl^ i(3?£$H. xi" 1 } - ^ g Tr - ifiitffi?} (A2) 



r 



Taking the variation respect to g™ k and k leads to the eigenvalue equations the matrix M. 
I 



sm ^ sin + + «* ^ C0S + <*-*) " ( cos ^ + cos ^ 



2N 2N 



N 



1 

2i 



N 



m.k 



(2j - 1j )g mk — 



4( ' Q Mg% k = 0, (A3) 



sin ^ sm J^(Cfe+i + Cn) + cos ^ cos ^(g™ k+1 + g^ k -i) ~ (cos ^ + cos ^)<# 



2N 2N 



N 



N 



±(£»-E«)3$_^M ff $ = 0, (A4) 



where £ = i8 s Gt/(4kGq). To solve the system of coupled equations it is convenient to write 

wb _ n >Jqk , < -iqk 

9 P ,k — 9 P e + g p e , 



(A5) 
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and similar equation for the m-component. The coeffi- 
cient q is to be determined from the boundary conditions 

n wb — n 

9p,0 — U : 



«a(0, x,N) = -(x + 2 cos((t>/N)) 



+ \/i(x + 2cos(^/7V)) 2 -l 



n wb — n 

From the boundary conditions we have — —gu and 
q = nl/N with I = 1, .., N — 1. We substitute the expres- 
sions ^3] in the equations for the eigenvalues and get 



4> b r nl / 

2 cos cos cos( — — cos — 

2N 2N y N' N 



9 d) b n 4> W o / 7rZ , „ . . , 



and the function 



2 s 2 



„2iV 



-2JV 



iV 



2 - s? iV _ sI 2N 



a 2N 



-2N 



N 



9 _ q 2N _ „-2JV 



■ln(«i(0,a: > JV)) 



5 

— \n(s 2 ((t),x,N)) 



The quantity e being (S b - E" , )/2^. 

Here we report the correction to the action when de- 
coherence and spin orbit are present as generalization of We can finall y ex P ress the correction to the action as 
Eq. [|3 , let us define the following functions 



S!{<f>,x,N) - \ (2 + t^tttt 

2 \ cos(0/iV) 



cos(0/7V) 



-S wl {4>) = A(4>, Ch,N)- 3A(<j>, Ch + ( so ,N) (A7) 
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